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I. 



CONTROL SYSTEM OVERVIEW 



This thesis is intended to act as the preliminary ground 
work for the attitude control system of the Naval 
Postgraduate School satellite ORION. The background for 
ORION is detailed in Reference 1. The orbital control 
system deals with the motion of the center of mass in an 
orbit around the earth, while the attitude control system 
deals with the angular motion of a satellite around its 
center of mass. The attitude control system must 
incorporate attitude sensors, control laws, actuators or 
thrusters, the effect of disturbing forces and the dynamics 
of the satellite. 

ORION is being designed to be launched via the extended 
Get-Away-Special (GAS) cannister aboard the Space 
Transportation System (Space Shuttle) . ORION is currently 
under development and many of the systems which will 
determine its final characteristics have not been fixed. 
For the purpose of this thesis, ORION is considered to be a 
cylindrical satellite of uniform mass with dimensions of 
0.889 meters height and 0.4826 meters diameter. It has a 
mass of 113 kilograms. Figure 1.1 shows a diagram of ORION. 

The design of a control system is highly dependent on 
the mission of the satellite. Reference 2 provides a 
section on the analysis of mission requirements and the 
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Figure 1.1 ORION Diagram 
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possible control systems which will fulfill mission 
requirements. Tables 1 and 2 show the trade-offs that need 
to be made when choosing a type of control system [Ref. 
2 : pp . IV-12 , IV-12 ] . 

ORION is being designed to act as a payload bus for a 
variety of missions. General requirements are that the 
satellite fly in orbits from 200 kilometers to 2000 
kilometers. Attitude sensor accuracies will not be less 
than 2.0 degrees. The life of the satellite will depend on 
a mission payload versus fuel trade-off. For missions that 
do not require a high degree of pointing accuracy, spin 
stabilization provides the simplest means of attitude 
control. A spinning satellite will be less susceptible to 
nutation if spun about a maximum moment of inertia and will 
be stable if energy can be dissipated. The equations 
governing satellite motion are developed in Chapter III. 

The purpose of the control system is to maintain the 
satellite at a specified orientation in space. If the 
satellite is perturbed from this orientation by external 
factors such as gravity gradient torques, aerodynamic drag 
or nutation then the control system must sense the change 
and apply the control torques necessary to reorient the 
satellite. It must also be able to reorient the satellite 
to accommodate any change in the mission requirements. 
Figure 1.2 shows a diagram of the basic control system. The 
feedback loop and error determination process may occur 
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Figure 1.2 Attitude Control Block Diagram 
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either on board the satellite or via a satellite to earth 
station telemetry link. On-board processing requires 
preprogramming the reference angles into the processor and 
is limited to those periods when the control system needs to 
act only as a regulator and any disturbances are small. A 
telemetry link is required to control the satellite during 
major orientation periods and flight maneuvers. It also 
acts as a back-up for the on board processor. 

There are three distinct phases of the attitude control 
problem. They are: 1) launch, consisting of activities 
required to take the satellite into a preliminary orbit? 2) 
acquisition, in which the satellite's preliminary 
orientation and maneuvers are accomplished; and 3) mission 
operations, in which the mission requirements are carried 
out [Ref. 3:p. 12]. 

A . LAUNCH 

Control during launch is limited to the launch vehicle 
and has little relationship to the attitude control of the 
satellite. Of primary concern is the detachment of the 
satellite from the launch vehicle in some preliminary orbit. 
The satellite will have the angular momentum of the launch 
vehicle plus some velocity imparted from the separation 
process. ORION will be carried into orbit by the Space 
Shuttle and launched from a GAS cannister. Reference 1 
details the cannister and the launch process. The launch 
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vector is along an outward radial from Earth's center at 4 
feet per second. 

B. ATTITUDE ACQUISITION 

The acquisition phase commences after ORION is detached 
from the Shuttle. During this phase the control problem 
consists of determining the attitude of the satellite, 
spinning the satellite up to the desired spin rate, 
reorienting the satellite to align the maneuvering thruster 
with the velocity vector required to place the satellite in 
its mission orbit and maintaining the correct attitude 
during the orbit transfer. Once in its final orbit, the 
control system must reorient the satellite to the desired 
mission attitude. 

NASA prohibits the firing of thrusters over 5 pounds 
sooner than 45 minutes after separation. Thrusters under 5 
pounds cannot be fired closer than 200 feet to the shuttle. 
At a launch velocity of 4 feet per second, this means the 
satellite will drift for 50 seconds before beginning spin up 
for attitude orientation. During this drift time the 

primary forces acting on the satellite are aerodynamic drag 
and gravity gradient torque. Additional torques will be 
present if the satellite contacts the cannister side during 
the launch process. The environmental torques are developed 
in Chapter III. 

The attitude of the satellite can be predicted during 
the drift period by applying these torques to the basic 
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satellite equations of motion. Some assumptions must be 
made concerning the initial conditions at the time of 
launch. These are: 

(1) The satellite does not impact the cannister during 
the launch process. If this happens the satellite 
will tumble and the attitude cannot be predicted. 
The motion would have to be observed by the Shuttle 
crew and reported. 

(2) The satellite is aligned as shown in Figure 1.3 with 
the z axis pointing towards the earth and the x axis 
in the direction of flight. The initial values of $, 
T and 0 are zero. 

(3) The aerodynamic drag force acts only against the 
satellite's center of pressure and opposite to the 
direction of flight. 

After the satellite is 200 feet from the shuttle, it is 
spun up to 60 rpm and the actual attitude determined by the 
horizon and sun sensors. A description of the horizon and 
sun sensors is given in Chapter II. Once the orientation of 
ORION is determined, control torques may be applied to 
orient it to any desired attitude. During this time system 
tests may be performed to verify the correct operation of 
the satellite. Details concerning control torques are found 
in Chapter III. 

The next portion of the acquisition phase is to orient 
ORION for orbital maneuvering. It is unlikely that the 
Shuttle will be able to launch ORION into its final orbit. 
For its first mission ORION will be required to extend four 
magnetometers from the body. This deployment will have two 
major effects. It will despin the satellite and it will 
increase the moments of inertia about the spin and 
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X \ 



transverse axes. Once the magnetometers are deployed, the 
satellite will need to be oriented to align the propulsion 
thruster with the velocity vector which will maneuver ORION 
into its final orbit. Details on the final orbit are not 
yet available. 

C. MISSION 

The third phase of the attitude problem is the control 
of the satellite while it is performing its mission. 
Generally this requires the control system to act as a 



regulator 


to maintain 


the correct 


orientation. 


If 


the 


required 


attitude is 


fixed, the 


control loop 


may 


be 



programmed into a micro-processor and function autonomously. 
The telemetry link would be used to monitor the processor 
and act as a back-up. Control by the earth station may be 
necessary to orient the satellite for maneuvers or to 
correct for a large disturbance. Figure 1.4 summarizes the 
control problem in algorithmic form. 

D. ASSUMPTIONS 

The purpose of this thesis is to form a model of ORION 
and simulate its motion during the three phases of the 
control problem. It will model the satellite's equations of 
motion, the environmental disturbance torques, and the 
control torques and study the interactions between these 
factors. As stated earlier, ORION is still in the early 
stages of design. Many of the systems which determine the 
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Figure 1.4 Attitude Control Algorithm 
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flight characteristics have not been developed. It is 
necessary therefore to make a number of assumptions which 
simplify the simulation problem. As the design of ORION 
matures assumptions may be removed and the model updated. 
This will help in determining design trade-offs and increase 
the accuracy of the model. 

The first major assumption is that the satellite has a 
uniformly distributed mass. This locates the center of mass 
at the origin of the body aligned axis as shown in Figure 
1.1. The moments of inertia are dependent on the deployment 
of the magnetometers, but in the launch configuration I a is 
3.239 kg-m 2 and I-j- is 8.916 kg-m 2 . As the mass distribution 
changes because of component location or changes in fuel 
weight, payload, etc. , these moments of inertia will change. 
If the moments of inertia are located along the body aligned 
axes, the cross products of inertia are zero. The second 
major assumption is that the orbit will be circular. This 
results in a constant angular velocity w 0 around the earth. 
The radius of orbit also remains constant throughout the 
orbit . 

The model of the satellite will first be constructed 
assuming the satellite is a rigid body with no internal 
energy effects such as sloshing fuel tanks. It will also be 
assumed that there will be no significant changes in the 
moments of inertia because of the control thruster burns. 
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Energy effects can be taken into account later in the 
simulation process. 

The model will also assume that the spin rate of the 
satellite will remain constant. In actuality the drag force 
on the satellite will cause the spin rate to decrease 
significantly in approximately 50 days [Ref. 4:p. 34]. This 
will be ignored in the model since if a significant 
reduction is noticed, then the spin thrusters can easily be 
fired to spin the satellite without affecting the 
orientation. 

When a satellite is in orbit, it is subjected to a 
variety of environmental forces and moments. These 
environmental disturbances include gravity gradient torque, 
aerodynamic drag, solar radiation pressure, magnetic field 
torques and micrometeor and radiation impacts. Gravity 
gradient and aerodynamic drag effects are dependent on the 
radius of the orbit. As the radius increases these effects 
are decreased. Solar pressure and magnetic torques are 
important over a long period of time, but can be ignored 
when looking at the short term motions. The primary effect 
of these environmental disturbances is noticed in the 
nutation and precession of the spinning satellite over time. 
This model will assume that the total torque on the 
satellite will be the sum of gravity gradient torque, 
aerodynamic drag torque and applied control torques. The 
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model will allow any additional torques to be applied at the 
simulator's discretion. 
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II. ATTITUDE SENSORS 



Before control can be applied to a satellite, the actual 
orientation in space must be determined. This is done by 
the attitude sensor package. Many devices and methods can 
be used as attitude sensors. These include rate gyroscopes, 
star sensors, sun sensors, earth limb or horizon sensors, 
and magnetometers. This chapter will describe how 
magnetometers, sun sensors and horizon sensors operate. 
These three types of sensors are currently planned for use 
in ORION. 

A. SUN SENSORS 

Sun sensors are used to measure the angle of the sun 
relative to the spin axis of the satellite. The sensor 
consists of a pattern of photoelectric cells and a housing 
that limits the angle at which the sun's rays contact the 
cells. Figure 2.1 shows a schematic of the sun sensor [Ref. 
3 :p. 163]. 

The field of view slit is oriented parallel to the x-z 
plane and perpendicular to the spin axis. The angle 8 de- 
termines the amount of sunlight on the photo cells which in 
turn produces a given amount of voltage. As the spin axis 
moves in inertial space, usually due to a torque, 8 changes 
and the voltage from the photo cells changes. The accuracy 
of the sun sensor is affected by the quality of solar cells 
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Source: [Ref. 3:p. 163] 

Figure 2.1 Sun Sensor Schematic 
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used, the presence of any highly reflective components of 
the satellite and the location of the sensor on the body. A 
typical accuracy is 0.01 degrees [Ref. 3:p. 17]. The rate 

of change of 3 is determined by calculating the time rate of 
change of the output voltage. 

Sun sensor data can be used to calculate nutation and to 
calculate the yaw angle ¥ . The method for determining 
nutation is described later in this chapter. ¥ is deter- 
mined by using information from the sun sensor and the 
horizon sensors. The motion of the spin axis is a combina- 
tion of both $ and T . If $ is known from the horizon 
sensors then for small angles, 

V 2 = g2 _ $ 2 (2.1) 

and 

• • • 

v = 3 - $ (2.2) 



B. HORIZON SENSORS 

Earth limb or horizon sensors are IR sensors with a set 
field of view that measures the difference in temperature 
between space and the earth's limb. Space is considered to 
have a temperature of 4 degrees Kelvin and the Earth has a 
temperature of 300 degrees Kelvin. As the sensor's field of 
view sweeps from space onto the earth a pulse is started. 
The pulse ends as the FOV goes from earth back onto space. 
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Two sensors are located equidistant from the center line 
of the satellite as shown in Figure 2.2. As the satellite 
spins, the two sensors create a pair of pulses as they cross 
the earth. The pulse width is determined by 

PW = R/2h0 (2.3) 

where R is the earth chord length (Figure 2.3), h is the 
altitude of the satellite and 0 is the spin rate of the 
satellite. 

• • • 
The sensors can determine $, $ , and 0. The spin rate© 

can be computed by measuring the time between successive 

pulse starts. $ is determined by comparing the pulse widths 

of sensor 1 and sensor 2 . When the pulse width of sensor 1 

is less than the pulse width of sensor 2, $ is positive. 

When the pulse width of sensor 1 is greater than that of 

sensor 2, $ is negative. This relationship is shown in 

Figure 2.3. The value of $ is computed by comparing the 

magnitudes of the pulse widths with those of $ = 0 and using 

the methods in Chapter 11 of Reference 5. $ is determined 

by comparing the pulse widths of successive pulses over 

time. 

The accuracy of the sensor is dependent on the quality 
of components, and the presence of infrared sources other 
than the earth. Masks are used to block the effects of the 
sun and avoid the saturation of the sensors. The satellite 
itself may act as a noise source, however shielding and 
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Figure 2. 




to. 



Y 



2 Horizon Sensor Location 



30 




SENSOR TRACK 
FOR $ > 0 





StNSOR TRACK 
FOR (j) = 0 



SENSOR TRACK 
FOR $ < 0 



Figure 2.3 Relationship Between Earth Chord and Roll Angle 
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knowing the infrared signature can reduce the effects. 
Pages 365-409 of Reference 5 give computational methods used 
to eliminate errors and determine the accuracy of the 
sensors. In general the accuracy of a horizon sensor can be 
assumed to be 0.1 degrees [Ref. 3:p. 17]. 

C. MAGNETOMETERS 

Considering the accuracy in attitude measurements 
required for ORION, the sun and horizon sensors will suffice 
for the sensor package. Since magnetometers are going to be 
used in the mission, they can be used to provide additional 
attitude inputs. These inputs can be compared with those 
obtained from the sun and horizon sensors to help eliminate 
noise and errors in measurements. 

For attitude sensing purposes, three mutually orthogonal 
magnetometers are required. These may either be internally 
mounted or deployed on booms [Ref. 3:p. 250]. The 
magnetometers consist of a coil around a ferromagnetic core 
of known permeability and area. When in the presence of a 
magnetic field, a voltage is produced that is proportional 
to the field component along the solenoid axis and is given 
by 



V = -ANy(dBj_/dt) (2.4) 

where is the field component [Ref. 3:p. 181]. By 

combining the field components from the three magnetometers, 
the solenoid axis directions can be computed. If the 
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solenoid axes are aligned with the body axes, the attitude 
angles can easily be found. 

Errors in the magnetometers may be the result of the 
model used to predict the earth's magnetic field, the 
influence of the satellite's electrical components, 
crosstalk between the magnetometers or a misalignment of the 
solenoid axis. At high orbital heights, the exact magnitude 
and direction of the earth's magnetic field is not well 
known and the model used to predict the magnitude and 
direction of the flux lines may include considerable errors. 
At altitudes above 1000 km, the field strength decreases as 
the inverse cube of the distance and the magnetic field of 
the satellite may dominate. Crosstalk involves the 
influence of one magnetometer on the others as a result of 
eddy currents that are produced. Misalignment may occur if 
the solenoid axis is not aligned with the body axis. This 
may occur upon deployment or if there is some deformation of 
the boom during a spinning maneuver. [Ref. 3:pp. 
181,182,250] 

D. NUTATION 

All sensor measurements described so far have assumed 
that the satellite was not nutating. If nutation occurs, 
the sensor data will exhibit an oscillatory behavior. The 
frequency of oscillation is the same as the body nutation 
rate u^. This motion can be used to calculate the nutation 
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angle directly from sensor data. The derivation and details 
are given on pages 539-548 of Reference 5. 

If the nutation is monitored by a sun sensor, then the 
maximum nutation angle is given by 

n = 6B/2R g (2.6) 

where 6g is the range of nominal sun vectors and R is the 
ratio of observed sun angle variations to the maximum 
nutation amplitude. R^ is dependent on the satellite's 
moments of inertia and the location of the sun sensor [Ref. 
3 : p. 250 ] . 

Knowing the nutation angle and the attitude angles, the 
satellite can determine when to fire control thrusters to 
reduce nutation and reorient itself. This can either be 
done autonomously or via the telemetry link. Examples of 
active nutation control and orientation are given in Chapter 
V. 
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III. ATTITUDE DYNAMICS AND CONTROL 



A. ATTITUDE DYNAMICS 

Newton and Euler first described the motion of a system 
of particles in inertial and rotational frames of reference. 
These are described and detailed in a number of 
publications. This chapter will present a summary of these 
derivations as they relate to modelling ORION'S motion. 

The satellite operates with respect to a reference frame 
which rotates around the earth in the orbital plane. As was 
shown in Figure 1.3, the x axis is in the direction of 
flight, the z axis points towards the earth and the y axis 
is perpendicular to the orbital plane. 

The orbital axes are related to the fixed reference 
frame of the earth by the direction cosine matrix C. C is a 
linear combination of Euler rotations defined as 



C = 




where F b is the body reference frame and F^ is the inertial 
reference frame [Ref. 5:p. 23]. The three rotations are 



C(<j)) 



10 0 

0 cos <J> sin 4> 

0 -sin <i> cos <P 



(3.1) 
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C(0) 



C(f) 



cos Q 0 -sin 0 

0 10 

sin 0 0 cos 0 _ 

cos ¥ sin H' 0 

-sin y cos Y 0 

0 0 1 

J 



(3.2) 



(3.3) 



The most commonly used sequence is the 3-2-1 sequence which 
results in 



C 



cosy cos 0 



sin^coscj) sin^sinif) 

-siny cos0sincj) +cos l l' sirf) coscj) 



-sin'i'cos 0 



cos’l'cosij) cos'Fsincj) 

+sinT sin0sincj> -cosj) sirtF sin0 



-sin0 



-sin<f>cos0 



COS0 COS(j) 



J 



(3.4) 



The satellite orientation is related to the orbital axes 
by the body fixed coordinate axes (x,y,z). The deviations 
of the body axes and the orbital axes are given by the 
angles (phi) , V (psi) , and 0 (theta) . 

There are six Euler equations that govern satellite 
motion [Ref. 5:p. 32]. They are: 



p 


= oo *p + f 


(3.5) 


p 


= mv 


(3.6) 


*c 


= V 


(3.7) 


h 


= -oo*h c + T 


(3.8) 


h 


= ICO 


(3.9) 


C 


= -co*C 


(3 . 10) 
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where oj * is defined as 







0 


-a) z 


CO 


( ii * 


= 


N 

3 


0 


-0) 






-U) y 


03 X 


0 



(3.11) 



Equations (3.5), (3.6), and (3.7) are the translation 

equations of motion and are not discussed in this thesis. 
Equations (3.8), (3.9), and (3.10) apply to the attitude 
control of ORION. 

The moment of inertia about an axis is the sum of the 
masses times the square of the distance of the mass from the 
axis. The moment of inertia matrix I is given as 







Ix 


Ixy 


Ixz 






i = 




lyx 


z y 


H 

N 




(3.12) 






Izx 


H 

N 


!z 






where I x = / x 2 dm, 


: 


■y = 


/y 2 dm, 


Iz ' 


= /z 2 


dm. Since ORION 


is a symmetrical 


body with the 


body 


axes 


centered at the 



center of mass, the cross-products of inertia, I X y, lyz an< ^ 
I zx are zero. 

Expanding Equation (3.8) into x, y, and z components 
yields 



h x = 


u z h y u y^z 




T x 


(3.13) 


hy = 


w x h x “ w z h x 


+ 


T y 


(3.14) 


h z = 






Tz 


(3.15) 
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Expanding Equation (3.9) and substituting into (3.13), 



(3.14), and (3.15) and solving for T yields 



x — Ix^x ( I z -I y) wyw z 



= I 



yO)y + (I x lz) U X W Z 



T z — Iz u z + ( ly” Ix) ^y^x 



(3.16) 

(3.17) 

(3.18) 



is the sum of the external torques on the satellite in 
the ith axis. These torques include environmental torques, 
control torques and any other torques which may affect the 
satellite. Torques will be discussed later in the chapter. 

Since ORION is symmetric about the x and z axes, I x = I z 
= 1^, the transverse moment of inertia. The y axis is also 
the spin axis and Iy = I a . The angular velocity vector a) 
decomposes into its Euler angle components: [Ref. 4:p. 30] 



03 . 



= 4> - Tco f 



W y — 0 Wq 



03 



= V +4)03, 



(3.19) 

(3.20) 

(3.21) 



ORION will spin such that the spin axis is in the negative y 
direction. This means that 



03y 



— 03o 0 



(3.22) 



In the absence of external torques the motion equations 
reduce to 



It w x — ( ^a *t) wy«z 

I t 03y = 0 



(3.23) 

(3.24) 
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(It- I a) 0J y aj x 



(3.25) 



It- co, 



t w z 



Equation (3.24) implies that Wy is a constant. Differen- 
tiating Equations (3.23) and (3.25) and substituting for w x 
and aj y yields 



■■ . ‘W 'W 2 

00 + ~ 00 00 =0 

x T 2 y x 



(3.26) 



" . ‘W ‘W 2 

03 + ~ 03 03 =0 

z 2 y z 



(3.27) 



For stability, the coefficients of 03 x and 03 z must be 
positive. The y axis is the maximum moment of inertia when 
I a > It and is the minimum moment of inertia when It > I a . 
If there is no energy dissipation, the satellite is stable 
if spun about either the maximum or minimum moment of 
inertia. If there are external torques or energy 

dissipation means, the satellite is stable only if spun 
about the maximum moment of inertia [Ref. 4:p. 112]. 

Since there are external torques present, these must be 
taken into account when calculating the attitude of the 
satellite. Substituting for co in Equations (3.16), (3.17), 

and (3.18), 

T x = I t (4>- V) + (Ifla) ( w o-0) (Y+o-b*) 

= I t 4 )_ It aj o^ + ( I t“ I a) ( oj 0 y+u^-OT - 0 ) o 04 )) 

.. 2 • • • 

= !•£<{>+ (It“^a) ( w o -w o® ) 4> + ( (It“^a) ( ^o”® ) “It w o) ^ (3.28) 
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(3.29) 



T z i t ( V/ ) + ( ) (w o” ©) (<f> — ) 

•• • • 2 • • • 

= It^ +I tWo4> + ( I a _I t) (oj o 4> - oJo'i'-0cJ) + coo0'l') 

= ( I a — 1-^ ) ( to o 0 ojo ) M / '*" ^*( 3-a ( too - 0) "*’-'-t (A) o) 4> 

= I t 'F+(It -I a) ( w o"w o 0)r( ( I t -I a) (too-0)- I t^o) 4> (3.30) 



Solving for the second derivatives of the Euler angles 
yields 

4> = (T x -(I t -I a ) (tOo-co o 0)4>-((I t -I a ) (co o -0)-I t to o )'i')/I t (3.31) 

Y = (\ - (It -I a) (uo““o©)^ + ((It- I a) (w o -0)-ItWo)i)/It (3.32) 

0 = -T y /I a (3.33) 

If <o 0 and 0 are assumed to be constant the equations can be 
written 

• • • 

cp = T x - A <f> ~ B (3.34) 

V = T z - Ay + B 4 , (3.35) 

0 = - Ty/I a (3.36) 

where : 

A = (1 ~ I a/ i t ) (too - 0)too 

B = (1 “ la/^t) (w o ”0) — w o 



The block diagram form of the equations is given in Figure 
3.1. As can be seen from Figure 3.1, <J> and are cross- 
coupled through their rotation rates. 
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Figure 3 . 1 System Block Diagram 




The total angular momentum vector H can be expressed as 



H - Ia^y + I t w t (3.37) 

2 2 2 . 

where wt = w x + w z an< * lies in the x-z plane. rotates 

relative to the satellite with the angular velocity A where 

A = (I a /It - l)wy (3.38) 

The body nutation frequency is defined as 



w n = I a/ I t 



(3.39) 



The motion of the satellite consists of the satellite 
rotating about its spin axis and the spin axis rotating 
about H with the angular velocity co n [Ref. 4:p. 115]. The 

nutation angle n i s defined as 



. -1 ,^t w t x 

H = sin (— — ) 



(3.40) 



Figure 3.2 shows the relationship between Wn , U y and H. 

The inherent stability of the satellite can be found by 
analyzing the transfer functions of the satellite. Setting 
T x and T z equal to zero and dividing by 1^. yields 

• • • • • 

(j)+(l-I a /I t ) (ojQ-ejajJ) +((1-I a /I t ) (u o -0)-o) o )’i' = 0 (3.41) 

V + (1-I a /I t ) (Uo-0)w^ -((1-I a /I t ) (w o -0) -u> 0 ) 4> = 0 (3.42) 

Noting that -A = (1-I a /It) (w o -0) and taking the Laplace 
transforms of the equations yields 



<fi(s)s 2 - A(jj 0 4>(s) - (A +co Q ) s 'i'(s) = 0 



(3.43) 
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a 



Figure 3.2 Nutation Cone 
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T(s)s 2 - Xo) 0 'F(s) + ( A+u)o) s<j)(S) = 0 (3.44) 

Solving for ^(s) and substituting 

2 2 

( A+io ) s 

<Hs) (s — Ato ) - ^ (J)(s) = 0 (3.45) 

° (s -Ao 3 q ) 

(J)(s) (s 2 - X(JJ Q ) 2 + ( A+Wq) 2 s 2 4>( s ) = 0 (3.46) 

(j,(s) (s 4 +( (A +03 o ) 2 -Aa3 0 ) s 2 + A 2 co q ) = 0 (3.47) 

Partitioning 



.2 2 
A o) Q 

cj)(s) = 1 + — 2 — j 



s (s +((A-Hjj ) -Aw )) 
o o 



The open loop transfer function is 



(3.48) 



A p (3.49) 

s 2 (s 2 + ( (A+a) Q ) 2 -Ato o ) 

The transfer function has two poles at the origin and two 
poles on the imaginary axis at +j J (A + w 0 ) 2 - Aoo 0 and 

-j y(A + “o) Aooo* Therefore, the stability of the satellite 

depends on spin rate, inertia ratio and orbit altitude. 
Table 3.1 shows the root locations for different 0 and ay. 

Because there is no cubed term in the transfer function, 
the roots will always be on the imaginary axis. This 

implies that the satellite will tend to reach a new steady- 

state value for <p and V if perturbed, rather than return to 
an initial state. 
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TABLE 3.1 



ROOTS OF THE TRANSFER FUNCTION 
(to 0 = -001184) 



CO y 


0 




-6.283 


0 .3635 


3 




1.4882 


-3 


-2 . 094 


0.3635 


1 




1.4882 


-1 


-1.047 


0.3635 






1.4882 


— 



Root Location 
A on Imaginary Axis 



999 


4 . 00 


067 


3.066 


333 


1.334 


022 


1.022 


6664 


.667 


5111 


.5106 



B . ENVIRONMENTAL TORQUES 

Before the motion of ORION can be accurately modelled, 
the possible external disturbances must be studied. The 
major disturbances experienced by a satellite in low earth 
orbit are aerodynamic drag and gravity gradient torques. 

Aerodynamic drag is a result of the motion of the 
satellite through the upper atmosphere. It can be expressed 
as 



D 




p V 2 C D A 



(3.50) 



where : 

p = density of air at orbit altitude 

V = velocity of satellite along orbit 

Cp = drag coefficient 

A = surface area over which the drag acts. 
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The drag on the body can be represented by a single force 
acting on the body's center of pressure. If this center of 
pressure is located away from the center of mass, a torque 
is developed. 

The velocity of the satellite is 



V = 



w 



w o R 



o — Pe/-^ 

V 2 = y e /R 3 *R 2 = i-ig/R 



(3.51) 

(3.52) 

(3.53) 



where : 

y e = gravitational constant 

R = distance of satellite from earth's center. 

The density of air decreases as the altitude increases. At 
a 200 km orbit the air density is approximately 7.0 x 10"^ 
kg/m 3 . The surface area is the cross-sectional area of the 
satellite; for ORION, A = 0.429 m 2 . Drag coefficients vary 
with shape, material smoothness, etc. A typical value for 
this type of satellite is Cp = 2. The drag can be expressed 
as 



D = 



14 

1.71 xlO P 
R 



(3.54) 



The center of pressure for ORION is assumed to be 
located 0.076 m along the positive y axis. Since the drag 
acts opposite to the direction of flight, the total drag 
torque on ORION can be expressed as 
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D * L * Cos Y 



(3.55) 



T D = 



where L is 0.076 m. 

The second major disturbance is due to the gravity 
gradient that exists as the distance from the earth 
increases. The gravity gradient torque is defined as 



3p 

G = — et / (r xr) (r*R) dm 

R 

o 

where : 



(3.56) 



R = position vector of the satellite with respect 
to the earth; and 

r = coordinate of an elemental mass from the 
center of mass. [Ref. 4:p. 131] 



Using the 3-2-1 Euler rotation between the body axes and the 
orbit fixed axis, the gravity gradient torques are 



G x = 


3 u o [(^t *a) s ^ n $ 


cos <p 


cos 2 © ] 


(3.57) 


G y = 

G z = 


0 

3u)o[(I t -I a )sinO 


cos 0 


sin <P ] 


(3.58) 

(3.59) 



Thus gravity gradient torques depend on orbital altitude, 
moments of inertia and satellite orientation. 

The equations describing aerodynamic drag and gravity 
gradient torques can now be used in simulating their effect 
on satellite motion. The analysis of this effect is 
contained in Chapter IV. 
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C. COMPUTER MODELS 



The computer model was written in the Dynamic Simulation 
Language (DSL) . The program listing is contained in the 
Appendix. The program was intended to be as general as 
possible to allow it to apply to a variety of parameter 
changes. Of special interest are 1^, I a , orbital altitude, 
air density, and spin rate. 

The three Euler equations are integrated once to obtain 
• • • 

<j> , y and Q . These are integrated again to obtain cf>, T , and 
Q . The initial values depend on the situation being 
analyzed. Once the attitude angles have been calculated, 
they can be used to calculate w t , to n , \ and the nutation 
angle. 

Also included is the equation governing magnetometer 



deployment . 


It 


is 


implemented 


in 


the 


form of a 


ramp 


function. 


The 


starting time 


is 


the 


time when 


the 


magnetometer 


is 


to 


be deployed. 




The 


ending time 


is 



calculated by dividing the desired boom length by 0.0169. 
If no boom length change is desired the ramp function can be 
commented out of the program and the variable Lj-, set to the 
desired boom length. An Ly of zero indicates that the 
magnetometers are in the storage position. 

Changes in spin rates can be accomplished in two ways. 
One is to change the initial condition of G in a PARAM 
statement. This is useful when evaluating the effect of 
configuration changes at different spin rates. The second 
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method simulates the use of the spin-up thrusters. Ty can 
be a step, ramp, or pulse input, or any combination that 
approximates the actual thruster firing profile. 

In addition to the variables describing moments of 
inertia, spin rate, and thruster torque, many other 
components of the model are left as variable quantities in 
the program. These can be changed by PARAM statements in 
multiple runs to provide the widest range of possible 
parameter studies. 
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IV. PARAMETER STUDIES AND ACQUISITION SIMULATION 



A. SIMULATION MODEL 

The design of any modern system involves the 
investigation and analysis of how the system responds to 
inputs with a given set of system parameters. These 
parameter studies are used in design and cost trade-offs 
which will eventually lead to the optimum system design. 
One way to run parameter studies is to use scale models and 
testing chambers. These are primarily used in the final 
design stages. In preliminary design where there are a lot 
of variables computer simulations are best performed. 

The computer model described in Chapter III was used to 
analyze some parameter changes that may be associated with 
ORION. It is expected that this model will be used to 
determine the optimum configuration in terms of moments and 
magnetometer deployment. This chapter analyzes the changes 
in inertia ratio and spin rates as a result of a typical 
mission. It also analyzes the environmental effects of 
gravity gradient and aerodynamic drag. 

B. PARAMETER STUDIES 

The first study was to analyze how a change in the 
moments of inertia would affect the spin rate and stability. 
The current mission requires the deployment of four 0.907 kg 
(2 lb) magnetometers on lightweight booms. The booms were 



50 



assumed to deploy at a constant rate of 0.0169 m/sec to a 
maximum of 3.048 m (10 ft). The mass of the booms was 
considered to be negligible and did not contribute to the 
inertia of the satellite. The magnetometers are located in 
the x-z plane of the satellite as shown in Figure 4.1. 

The effect of deploying the magnetometers is to increase 
both I a and If The new moments of inertia can be 
calculated by finding the moment of inertia for the 
satellite without the magnetometers and then adding the 
moment of inertia for each magnetometer. Thus 

I a = 0.5m s r 2 + 4m m (l b + r) 2 (4.1) 

I t = 0. 083m s (3r 2 +h 2 ) + 2m m (l b +r) 2 (4.2) 

where : 

m s = mass of satellite minus the mass of 4 
magnetometers 

m m = mass of one magnetometer 
r = radius of the satellite 

h = height of the satellite 

l b = distance of the magnetometer from the outside 
of the satellite. 

Substituting in the values for ORION, the moments of inertia 
can be found at any boom length by 

I a = 3.203 + 3 . 628 ( l b + 0.242) 2 (4.3) 

I t = 8.81 + 1 . 8 1 ( l b + 0 . 242 ) 2 (4.4) 
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Figure 4 . 1 Magnetometer Deployment Diagram 
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Figure 4.2 shows this relationship between the boom length 
and the moments of inertia. I a becomes larger than I t at a 
boom length of 1.53 m. If energy could be dissipated, then 
with I a larger than I-j- the satellite would become stable 
about the spin axis [Ref. 4:p. 117]. 

Deployment of magnetometers also has an effect on the 
spin rate and could be used to change it without using the 
spin-up thrusters. Using the principle of the conservation 
of angular momentum, if co^ is assumed to be zero then 

H = ^a^y (4.5) 

and as I a increases tOy must decrease. Figure 4.3 shows how 
the spin rate changes with boom length. 

As noted earlier, the magnetometer deployment can be 
used to change the spin rate of the satellite. If the sat- 
ellite is spinning at 60 rpm (Wy = -6.283 rad/sec), a boom 
length of 1.097 m will slow the spin to 20 rpm (Wy = -2.094 
rad/sec). A 1.87 m deployment will slow the spin to 10 rpm 
(coy = -1.047 rad/sec). Adjusting the boom length to reach a 
specific spin rate is applicable if the mission does not re- 
quire the magnetometers a specific distance from the satel- 
lite. It is more probable that the boom must be a specific 
length and the satellite must spin at a specific rate. 

There are two ways to achieve this configuration. The 
first method is to deploy the magnetometer to the desired 
length, then adjust the spin as necessary with the 
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Figure 4 . 2 Change in Moments of Inertia with Boom Length 
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Figure 4.3 Change in Spin Rate with Magnetometer Deployment 
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thrusters. A second method is to despin the satellite with 
the thrusters then deploy the magnetometers and readjust the 
spin with the thrusters. If the change in spin and boom 
length is small, it may be possible to despin the satellite 
to a rate such that the subsequent deployment of the booms 
will despin the satellite to the required rate. 

Since the firing interval of the thrusters is 
proportional to the amount of fuel used, it can be used to 
compare the two methods. To illustrate this, assume the 
satellite is initially spinning at 60 rpm with I a = 3.239 
kg-m 2 . A boom length of 2 m and a spin rate of 20 rpm is 
required. At 2 m, I a is 21.447 kg-m 2 . Using the first 
method, after deployment the spin rate is 

I . w , 
al yl 

I o 
a2 

(3.239) (-6.283) 

21.477 

-0.9476 rad/sec (4*6) 




Assuming the total torque of the thrusters is FxR = 0.2168 
n-m, the firing interval dt is 



dt = 



I a2 (2 - 094 -“y2> 
F xR 



(21.447) (2.094-0.9476) 
0.2168 

= 113.3 sec 



(4.7) 



56 



Using the second method, the firing interval required to 
despin the satellite to 20 rpm is 

(3.239) (6.283-2.093) 
dt = (0.216 51 

= 62.6 sec (4.8) 

Deploying the magnetometer to 2 meters causes the spin rate 
to decrease to 0.3152 rad/sec. The firing interval required 
to readjust the spin is 17 6.1 sec. The total firing 
interval for the second method is 238.7 sec, more than twice 
that of the first method. 

This indicates that the first method is more fuel 
efficient in deploying the magnetometers and achieving a 
specific spin rate. A disadvantage of both these methods is 
that if the booms are too flexible, adjusting the spin may 
cause the magnetometers to lag the boom-satellite attachment 
points and cause too great a stress in the booms. 

1 . Stability 

Whether a is changed by deploying magnetometers or 
by changing the internal mass distribution of the satellite, 
it will have an effect on the stability of the satellite. 
Figure 4.4 shows how a change in c affects X and w n . Since 
the roots of the open loop transfer function depend on X , 
they also change with a as was shown in Table 3.1. Figure 
4.5 shows the responses of the satellite to a 1.0 n-m pulse 
of 0.01 sec as a changes. This was done at three different 
spin rates of 60 rpm, 20 rpm, and 10 rpm. The advantage of 
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Figure 4.4 The Effects of Inertia Ratio on 
Nutation Frequencies 
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Figure 4.5 Satellite Response to a Short Duration Torque 
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increasing o is to make the satellite less susceptible to 
nutation. Figure 4.6 shows that the higher the spin rate 
and the larger the inertia ratio, the smaller the nutation 
will be for a given input torque. 

C. ENVIRONMENTAL EFFECTS 

The aerodynamic drag and gravity gradient torques were 
modelled using the equations in Chapter III. These 
equations were incorporated into the simulation to see what 
effect they would have on the satellite. 

Since the drag force is dependent on the density of the 
atmosphere, it decreases as the orbit height increases. 
This is shown in Table 4.1. At an altitude of 200 km the 
drag torque is 1.383 xlO -3 n-m and decreases to 1.762 x 10" 9 
at an altitude of 1000 km. 

TABLE 4.1 

TYPICAL VALUES FOR AERODYNAMIC DRAG 



ALT (km) 


Density (kg/m 3 ) 


Drag (n) 


T D (n-m) 


100 


6 xlO -7 


15.84 


1.204 


200 


7 xlO -10 


. 0182 


1.383 xlO -3 


400 


1 x 10 -11 


2.523 x 10" 4 


1.917 x 10 -5 


1000 


1 x 10" 15 


2.318 x 10 -8 


1.762 x 10 -9 



Gravity gradient torque is also dependent on orbital 
altitude since it depends on w 0 . Unlike drag, however, 
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Figure 4.6 Nutation Angles for Different Spins and Inertias 
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gravity gradient torque is not a constant but depends on the 
moments of inertia and the orientation of the satellite. At 
an altitude of 200 km the maximum torque is 3.99 xio -9 n-m. 
This decreases to 3.72 x 10~ 15 n-m at 1000 km. 

An orbital altitude of 200 km was used to determine the 
effect environmental forces have on ORION. Inertia ratios 
of 0.364 and 1.488 and spin rates of 10, 20 and 60 rpm were 
used to give the widest ranges of results. These results 
are shown in Figures 4.7, 4.8 and 4.9. At 10 rpm the 
maximum nutation is 2.053 x io -3 radians for a a of 0.364 and 
5.33 x 10~ 5 radians for a a of 1.488. At 20 rpm the nutation 
angles are 5.173 x 10 -4 radians and 1.17 x io -5 radians and 
for 60 rpm they are 6.0 x io“^ radians and 1.33 xl0“® 
radians. 

As can be seen in the figures, the nutation is 
oscillatory. This is because of the oscillatory nature of 
the environmental torques, principally the gravity gradient. 
When the satellite nutates to its maximum angle, the torque 
effect is reversed causing the nutation to decrease. When 
the nutation is zero, the torque effect is again reversed 
and the process repeats itself. The frequency of the effect 
is related to the spin rate of the satellite. This is seen 
in the 0 term of the torque equations. 

If the accuracy requirement of the attitude control 
system is small enough, environmental effects must be taken 
into account in the simulations. They will be ignored in 
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Figure 4.7 Nutation Due to Environmental Torques; 10 Rpm 
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Figure 4.8 Nutation Due to Environmental Torques; 20 Rpm 
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Figure 4.9 Nutation Due to Environmental Torques; 60 Rpm 
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further simulations described in this thesis since their 
effect is less than the 2 degree accuracy requirement of 
ORION. This will simplify the model without a large error 
in accuracy. 

D. ACQUISITION PHASE PROFILE 

After analyzing the effects of o changes and 
environmental torques, the next step was to simulate the 
motion of the satellite from launch to spin-up. Drag and 
gravity gradients effects were included in this simulation. 

The flight profile consisted of launching the satellite 
with no spin. All initial values of cf> , 'i and 0 were assumed 
to be zero. At 50 seconds a step input of Ty = 0.2168 n-m 
was applied. This caused the value of «y to increase. At 
144 seconds the input was removed. The value of coy was 
-6.29 rad/sec (60 rpm) . Figure 4.10 shows how a and Y 
change during this portion of the mission. 

During the drift period environmental torques cause the 
satellite to rotate 0.013 degrees in roll and 11.04 degrees 
in yaw. As the spin increases the a and 'H begin an 
oscillatory motion. The amplitude of the motion depends on 
the spin rate. As the spin rate increases, the satellite 
becomes more resistive to torque effects. When the spin 
rate becomes a constant, the amplitude of oscillation due to 
environmental torques becomes a constant although the scale 
in Figure 4.10 is too large to show this. At 144 seconds, <f) 
is 3.89 degrees and ¥ is 12.54 degrees. The nutation 
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Figure 4.10 Acquisition Phase Roll and Yaw 
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caused by environmental torques is 0.0015 degrees, much less 
than the acceptable nutation of 2 degrees. After 400 
seconds with no controlling inputs, <J> is 7.4 degrees, V is 
10.8 degrees and the nutation is 0.0015 degrees. 

From this point, the magnetometers could be deployed or 
control torques could be applied to reorient the satellite. 
Chapter V discusses the methods of controlling nutation and 
reorienting the satellite. 



68 



V. MISSION PHASE SIMULATION 



A. SATELLITE CONTROL 

The process of maintaining the correct orientation of a 
satellite consists of eliminating the precession and 
nutation effects caused by external or internal 
disturbances. The simplest case is when the spin axis and 
angular momentum vector are aligned and are perpendicular to 
the orbital plane. 

This presents a two part control problem. The first 
part is to reduce nutation by forcing equal to zero. The 
second part is maintaining the angle between the satellite's 
angular momentum vector H and the orbital angular momentum 
vector H 0 . The relationship between H and H 0 is 

H -H 

0 

Cos a - 

H | |H 

0 1 1 

The angle a can be obtained by precessing the satellite with 

control torques and then eliminating any resulting co t . 

to £ can be forced to zero by applying a control torque of 

magnitude equal to I-t^t i- n the direction of I-t w f w t can te 

• • 

calculated from measurements of <J>, ¥ , <f> and ¥ . The 

magnitude of torque required is determined by 

F * r * dt = It w t 



where: 
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F 



the impulsive force of the thruster 



r 



distance of the thruster from the center 
of mass 



dt 



time over which the thruster is fired 



An ideal thruster would provide an impulsive force at the 
correct time in the spin cycle and dt would be zero. The 
firing profile of an actual thruster can be approximated by 
the trapezoid model shown in Figure 5.1 [Ref. 3:p. 21]. t-^ 
is the time the thruster is commanded on, t 2 is the time the 
thruster reaches maximum force, t 4 is the time the thruster 
is commanded off and t 5 is the time the thruster force is 
zero. The total time dt is from t^ to t 5 . t c is the time 
at which the effective torque is applied. 

t 2 and t 5 can be determined from experimental or the 
manufacturer's data. Knowing t c , t 2 and t 5 , and t 4 can 
be determined by 



Tan At c = b/a 



(5.1) 



a cos At + b sin At 
c c 



(5.2) 



r 



where 



cos At 2 - cos At^ 



cos At,. - cos At^ 



a 





sin At 2 - sin At^ 



sin At,. - sin At 4 



b 
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Figure 5.1 Typical Thruster Profile (Trapezoid Model) 
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Figure 5.2 shows the relationships between the firing 
profile and timing. 

An even simpler approximation is the pulse function. 
This would occur if t 2 = ti and t 4 = t 5 . Knowing that dt = 
t 4 -t 2 , t c = (l/2)dt. Therefore 

t 2 = t c - ( 1/ 2 ) dt (5.3) 

t 4 = t c + ( 1/2 ) dt (5.4) 

Once the magnitude of control torque and the firing arc 
are known, the problem is when during the spin cycle to fire 
the thruster to eliminate The acceleration of any point 

(x,y,z) on the satellite is 

2 

a x = -X(o s + to n yA sin At (5.5) 

a y = (ojs-A)A(z cos At + x sin At) (5.6) 

2 

a z = -Zo) s + u) n y Acos (5.7) 

where A is proportional to the nutation angle [Ref. 4:p. 
116]. If an accelerometer is placed in the x-z plane, the 
axial acceleration reaches a positive peak when points 
towards the accelerometer. Placing a thruster 90 degrees 
from the accelerometer produces a torque opposite to “f If 
the accelerometer is located on the z axis, 

a y = (cog-AjAz cos At (5.8) 

and the peak acceleration is reached when day/dt = 0. 
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Figure 5.2 Relationship Between Spin Rate and 
Firing Times 
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Since 



then 



day/dt = -(oj s -A)AzA sin At 



(5.9) 



day/dt = 0 when t = it /A 

The value of acceleration at the firing point is 

a ymax = — (w s — A)zA (5.10) 

Using the value of t = tt/A as t c , the on and off times for 
the thruster can be calculated. If the spin rate is too 
fast, dt may encompass several revolutions of the satellite. 
To provide a torque in a given direction, however, the 
firing must be limited to an arc centered about t c . 
Therefore, the total firing time required (dt) must be 
divided into segments occurring at integer multiples of t c . 

For example, assume the nutation angle is 1.0 deg 
(0.0174 rad), I t is 8.91 kg-m 2 , W y is -6.283 rad/sec, and H 
is 20.35 kg-m 2 /sec 2 . Therefore 

I t oj t = H sin n 

= 0.355 kg-m 2 /sec 2 (5.11) 



Assuming F is 0.448 n and R is 0.242 m 



dt = 



Fr 



= 3.28 sec 



(5.12) 
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At 60 rpm the satellite makes one revolution per second. If 
the firing arc is 90 degrees, the firing time per arc is 

TT 

2 ~j— = 0.25 sec (5.13) 

Therefore, 13 firings will be required to eliminate the 
nutation. If the spin rate was 10 rpm, only 2 firings would 
be required. Each firing must occur at ay max to avoid 
adding energy to the system. 

This method of active nutation control was simulated 
using the computer model . Nutation was started by inputting 
a 5 n-m pulse of 0.01 seconds. After 0.03 seconds the 
control torque could be applied. The minimum required 
nutation was set at 0.001 radians to avoid computational 
errors. The firing arc was limited to 45 degrees. 

To apply the correct control torque, the program 
calculated the value of ay which would cause the 
acceleration peak to be in the center of the firing 
interval. When the calculated value of ay was greater than 
this value, T z was set equal to 2Fr/dt where dt is the 
firing arc time period. As long as the nutation was greater 
than 0.001, T z was applied at the peak positive acceleration 
point . 

This simulation was run for 3 different spin rates of 
10, 20 and 60 rpm at an inertia ratio of 0.363. It was also 
run for two different inertia ratios and a constant spin 
rate. These results are shown in 3 series of 3 figures. 
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Each series shows the nutation angle and applied control 
torque, the $ and 'F angles and the firing diagrams for a 
different spin rate. 

In Figure 5.3, the spin rate was 10 rpm. The initial 
nutation was 1.45 x 10 -3 radians. It was reduced to 0.001 
radians at 4.4 seconds. T z was applied to reduce the 

nutation to 0.001 radians. From Figure 5.4, <J> and y were 

reduced to steady-state values of 1.67 X 10 -2 radians and 
-1.66 x 10” 4 radians. The values of <f> and V oscillated about 
these values because the nutation was not completely 
eliminated. Figure 5.5 shows the firing diagram for this 
case. T z was 8.33 x 10” 2 n-m over a firing time of .6 sec. 

In Figure 5.6, the spin rate was changed to 20 rpm which 
resulted in an initial nutation of 7.28 x 10~ 3 . The first of 
3 control firings occurred at 2.2 seconds. Nutation was 
reduced after 12 seconds. In Figure 5.7, the values of <j> 
and y are 7.91 x io -3 radians and -2.1 xio -4 radians. From 
the firing diagram of Figure 5.8, T z was 4.1 x 10 -2 n-m over 
a pulse of .36 sec. Three pulses were needed for a total 
firing time of 1.08 seconds. 



In 


Figure 5.9, 


a spin rate of 


60 rpm was 


used . 


The 


initial 


nutation 


was 2 . 43 x 10” 3 


radians , 


the 


first 



application of control occurred at .7 seconds and the 
nutation was reduced after 27.5 seconds. $ and 'F from 
Figure 5.10 were 2.27 x io -3 and -2.0 x io"^ radians. From 
the firing diagram of Figure 5.11, 18 applications of T z 



76 



[ouut.ou m 




Figure 5.3 Nutation Control? ujy = 10 Rpm 
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Figure 5.4 Roll and Yaw Angles; = io Rpm 
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Figure 5.5 Firing Diagram; oiy = io Rpm 
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Figure 5.6 Nutation Control; = 20 Rpm 
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Figure 5.8 Firing Diagram; uy = 20 Rpm 
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Figure 5.9 Nutation Control; oiy = 60 Rpm 
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Figure 5.10 Roll and Yaw Angles; toy = 60 Rpm 
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were required. T z was 1.3 xio -2 n-m with a pulse width of 
.23 seconds. A total firing time of 4.2 seconds was 
required . 

Figure 5.12 shows the difference in control for the same 
spin rate of 20 rpm but for two different inertia ratios. 
When o = 0.363 it required 3 applications to reduce the 
nutation. When was increased to 1.12 the initial amount 
of nutation was less and only one application of control 
torque was required. 

As a result of these simulations, it can be seen that 
the slower spin rates require less firing times and number 
of pulses. This is because the firing arc time is larger. 
A low number of pulses also reduces the probability of a 
timing error which may cause nutation instead of eliminate 
it. In every case shown, <p reached a steady state value 
while y approached zero. This is a result of using T z for 
the input and control torques. 

Fewer pulses resulted in shorter firing times and 
quicker responses. This equates to a fuel savings. The 
model can also be run with different size thrusters by 
changing the value of F in the program. Larger firing arcs 
would also decrease the response times especially at higher 
spin rates. 

B. SATELLITE REORIENTATION 

The control system may be used to reorient the satellite 
as well as reduce the nutation. Reorientation can 
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Figure 5.12 Nutation Variation with^ ; coy = 20 Rpm 
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be accomplished either as part of the nutation reduction or 
as a separate control sequence. It was shown in Figures 
5.4, 5.7 and 5.10 that although the nutation was reduced, <P 

and 'F were not returned to their initial values. 

Reorientation during nutation control is a two-step 
process. The first step involves reducing the nutation 
angle by half its original amount. This will cause a new 
nutation cone to be established. When the spin axis moves 
180° along the new cone, control torques are applied to 
reduce the nutation to the minimum level. 

Because of the differences between nutation frequency, 
spin rate and A , the thrusters may not be aligned with the 
vector when the satellite reaches the 180° point on the 
nutation cone. The time required for the spin axis to 
travel 180° is 

At = n/wn (5.14) 

The satellite will rotate about its axis through an angle 

a = w y At (5.15) 

in this time period. The relative rotation of the vector 
is 



8 = A At (5.16) 

over the same time period. The two angles o. and 8 are equal 
only if a = 1.0. 
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Figure 5.13 shows the relationship between Wy, w n and X. 



Since 



X 



(a-1) Wy 

(I a /I t ~ l)0)y 



(5.17) 



and 



w 



n 




(5 . 18) 



The two angles a and 3 can be calculated as 



3 



tt (a-1) 

<i> 



(5.19) 



a 



it / a 



(5.20) 



Since over the required time period the thrust vector may 
not be aligned with the vector, the application of a 
control torque would add to the nutation. 

To overcome this difficulty, four sets of thrusters can 
be used. If the thrusters are located on the x and z axes, 
a combination of thrusters can be fired to make the 
resulting torque coll inear with the vector. 

This is accomplished in the program by saving the time 
at which the nutation is halved as FTIME. When the program 
time reaches FTIME plus At, 3 is calculated and the correct 
combination of T x and T z used to reduce the nutation to an 
acceptable level and restore <J) and ¥ to their initial 
values . 
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The second method is to reorient the satellite assuming 
there is no nutation. The degree of reorientation was given 
earlier in the chapter as 

H -H 

cos a = ThJThT 

Knowing the angle a and H 0 , the H required to make a zero 
can be calculated. This gives the required A control 
torque equal to Itw ^ is applied to start the satellite 
nutating. When A t is equal to iT/w n a control torque equal 
and opposite to It^t is applied to reduce the nutation. 
This is done in a manner similar to that of reorienting the 
satellite with nutation. 
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VI. CONCLUSIONS 



This thesis derived the equations which describe the 
rotational motion of a satellite. It also derived the 
equations describing the effects of gravity gradient torque 
and aerodynamic drag. These equations were then used to 

construct a computer simulation model . The model was used 
to analyze how changes in moments of inertia and spin rates 
affect the motion of the satellite. 

It was shown that extending magnetometers on booms would 
cause I a and 1^. to increase. This changed the inertia ratio 
a. The spin rate of the satellite could then be changed by 
changing the boom length of the magnetometer. As c 
increased, the stability of the satellite was changed. 

A high value of a made the satellite resistive to nutation 
by external forces. 

The effects of gravity gradient and aerodynamic drag 
were analyzed. The amount of nutation caused by environment 
depended on the spin rate and a. The maximum amount of 
nutation occurs at a low spin rate and low a . The nutation 
is on the order of 2.053 x 10 -3 radians for a oof 0.363 and 
a 10 rpm spin rate. The effects of the environment need 
only to be taken into account in the motion analysis if the 
accuracy required by the control system is small. The 
accuracy for ORION is to be within 2 degrees so gravity 



92 



gradient and drag were only applied during the example 
flight profile. 

The control system that is used is an active nutation 
control. An example of how the control system would be used 
was given. Assuming a firing arc of 45°, 4.4 seconds were 
required to reduce the nutation at a spin rate of 10 rpm, 12 
seconds for a spin rate of 20 rpm and 27.5 seconds for a 
spin rate of 60 rpm. The procedure for reorienting the 
satellite to any desired attitude was also described. 

Although the computer model was written to simulate the 
satellite as closely as possible, many areas need more 
detailed modelling. The satellite was assumed to be a 
cylinder of uniform mass. The actual mass distribution 
needs to be determined and the moments of inertia 
calculated. This would change the inertia ratio and 
possibly cause some cross-products of inertia to develop. 

Another factor that was not included in this model is 
the effect of liquids slosh in the fuel tanks. This could 
cause internal torques that adversely affect the motion. 
These tanks could also be designed to dissipate energy. If 
this could be done, then when o was greater than one, the 
satellite could be made stable. 

The firing profile of the control thruster was assumed 
in the model to be a pulse function. Entering the profile 
as a trapezoid would increase the accuracy of the model. 



93 



When real thruster data becomes available it could then be 



entered. 

As has been stated before, many of the variables 
program can easily be changed in parameter runs, 
allows a great deal of flexibility in analyzing 
changes as they occur. 



in the 
This 
design 
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APPENDIX 

ATTITUDE ANGLE CONTROL PROGRAM 



TITLE ATTITUDE ANGLE CONTROL PROGRAM 

* THIS PROGRAM PROVIDES THE SIMULATION MODEL FOR THE CONTROL OF THE 

* ORION SATELLITE. IT USES THE EOUATIONS DERIVED IN CHAPTER 3 TO MODEL 

* SATELLITE MOTION AMD ENVIRONMENTAL TORQUES. IT CALCULATES THE ANGULAR 

* FREQUENCIES AND THE ATTITUDE ANGLES. IT CALCULATES THE MAGNITUDE OF 

* THE ENVIRONMENTAL TOROUES BASED ON ORIBTAL PARAMETERS. IT ALSO 

* COMPUTES THE NUTATION; ANGULAR MOMENTUM VECTOR AND THE CONTROL TORQUES 

* NECESSARY FOR THE ACTIVE NUTATION CONTROL SYSTEM. CHANGES TO THE 

* PROGRAM CAN BE MADE VIA PARAM STATEMENTS OR DIRECTLY IN THE LISTING. 



* VARIABLE DESCRIPTION 



x t, 

* F 

* R 

* IT 

* IA 

* ALT 

* RO 

* ROE 

* WO 

* D 

* wx 

* wz 

x WY 

* FHID 

* PSID 

* THD 

* PHI 

* PSI 

* TH 

* TDX 

* TDZ 

* TD 

* GX 

* GZ 

* GG 

* IX 

* TZ 

X TY 

* WT 

* HT 

* HS 

* H 

* NU 

* SIGMA 

* LAND A 

* WN 

* AY 

* FARC 

* TC 

* AON 

* FTIME 



DISTANCE FROM THE CENTER OF PRESSURE TO THE CENTER OF MASS 

THRUST FORCE OF ONE CONTROL THRUSTER 

RADIUS OF THE SATELLITE 

TRANSVERSE MOMENT OF INERTIA 

MOMENT OF INERTIA ALONG THE SPIN AXIS 

ORBITAL ALTITUDE OF THE SATELLITE 

RADIUS OF THE EARTH 

DENSITY OF AIR AT THE ORBITAL ALTITUDE 
ANGULAR VELOCITY OF THE SATELLITE WRT EARTH 
DRAG FORCE 

ANGULAR VELOCITY ABOUT THE X AXIS WRT EARTH 

ANGULAR VELOCITY ABOUT THE Z AXIS WRT EARTH 

ANGULAR VELOCITY ABOUT THE Y AXIS WRT EARTH 

SATELLITE ROTATION ABOUT THE X AXIS 
SATELLITE ROTATION ABOUT THE Z AXIS 
SATELLITE ROTATION ABOUT THE Y AXIS 
ROLL ANGLE 
YAW ANGLE 
PITCH ANGLE 

DRAG TORQUE IN THE X DIRECTION 
DRAG TORQUE IN THE Z DIRECTION 
TOTAL DRAG TOROUE 

GRAVITY GRADIENT IN THE X DIRECTION 
GRAVITY GRADIENT IN THE Z DIRECTION 
TOTAL GRAVITY GRADIENT TORQUE 
CONTROL TOROUE IN THE X DIRECTION 
CONTROL TORQUE IN THE Z DIRECTION 
CONTROL TORQUE IN THE Y DIRECTION 
TRANSVERSE ANGULAR VELOCITY 

TRANSVERSE COMPONENT OF THE ANGULAR MOMENTUM VECTOR 
AXIAL COMPONENT OF THE ANGULAR MOMENTUM VECTOR 
MAGNITUDE OF THE ANGULAR MOMENTUM VECTOR 
NUTATION ANGLE 
INERTIA RATIO 

RELATIVE VELOCITY OF WT WRT THE SATELLITE ROTATION 

NUTATION FREQUENCY WRT INERTIAL SPACE 

AXIAL ACCELERATION OF THE SATELLITE DURING NUTATION 

FIRING ARC OF THE CONTROL THRUSTERS 

TIME OF PEAK POSITIVE AXIAL ACCELERATION 

ACCELERATION AT THE THRUSTER TURN ON TIME 

TIME OF THE THRUSTER FIRING 



* CANCEL = TIME FOR THE SPIN AXIS TO TRAVEL PI RADIANS 



CONST IT=8 .911 , IA=3 . 239 , R=0.242 ( F=0.448, NUREF=0 .001 ,L=0 . 076 
METHOD RKSFX 
FIXED FLAG 
INTEGER 

* INITIAL PROGRAM STATEMENTS: ORBIT HEIGHT, DENSITY, ORBIT VELOCITY 
INIT 
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ALT = 200.00 

ROE = 7.0D-10 

RO = (6378 . 0+ALT)*l . 0D03 

WO = SQR7(3.987D14/RO**3) 

D = 1.7D14*ROE/RO 
FLAG = 0 

INCON WXO=0 . 0 , WZO=0 . 0 , WYO=- 2 . 09 3 , PHIO=0 . 0 , PS 10=0 . 0 , THO=0 . 0 
DERIV 

WXD = (TX+TDX+GX- (IT-IA)*WY*WZ) /IT 

WZD = (TZ+TDZ+GZ- ( I A- IT ) *WY*WX) /IT 

WYD = TY/IA 

PHID = WX + PSI*WO 

PSID = WZ-PHI*WO 

THD = WY+WO 

WX = I MTGRL ( WXO , WXD ) 

WZ = INTGRL ( WZO , WZD ) 

WY = INTGRL (WYO, WYD) 

PHI = INTGRL (PHIO, PHID) 

PSI = INTGRL (PS 10, PS ID) 

TH = INTGRL (THO, THD) 

DYNAMIC 

* ENVIRONMENTAL TOROUE CALCULATION; USE WHEN DISTRUBANCE TORQUES ARE 

* REQUIRED 

* TDX = -D*L*SIN(TH) *COS (PSI ) 

* TDZ = D*L*COS(TH)*COS(PSI) 

* TD = (TDX**2+TDZ**2)**0.5 

* GX = 3.0*(WO**2)*(IT-IA)*SIN(PHI)*COS(PHI)*(COS(TH) )**2 

* GZ = 3.0*(WO**2 WlT-IA)*SIN(TH)*COS(TH)*SIN(PHI) 

* GG = (GX**2+GZ**2)**0.5 

■k-k-k'k'k-r.'k 

* USE THE NEXT 4 STATEMENTS WHEN NO DISTURBANCE TORQUES ARE REQUIRED 

TDX =0.0 
TDZ =0.0 
GX = 0.0 
GZ =0.0 

TX = 0.0 
TY =0.0 

* 

* STEP FUNCTION USED TO SET AN INITIAL NUTATION VALUE 

TZ = 5.0*STEP(0.01)-5.0*STEP(0.02) 

X 

* CALCULATION OF ANGULAR VELOCITIES, MOMENTUM VECTORS AND FREQUENCIES 

WT = (WX**2 + WZ**2)**0.5 
HT = IT*WT 
HS = IA*WY 

H = (HT**2 + HS**2)**0.5 
NU = ARSIN (HT/H) 

SIGMA = I A/ IT 

LAMDA = (SIGMA-1 . 0)*WY 

WN = SIGMA *WY 

X 

* CALCULATION OF THE FIRING TIME AND FIRING ARC 

AY= (WY- LAMDA) * R * COS (LAMDA*TIME) 

FARC = 0 . 392699/ABS (WY) 

TC = 3 . 14159/LAMDA 

AON = (WY-LAMDA) *R*COS (LAMDA* (TC- FARC) ) 

* THIS STATEMENT ALLOWS THE NUTATION TO BE PRESENT BEFORE CONTROL OCCURS 

IF (TIME. LE. 0.2) GO TO 15 

X 

* CONTROL CALCULATION AND FIRING. IF FLAG IS 0 THE CONTROL SYSTEM WILL 

* REDUCE THE INITIAL NUTATION BY HALF THEN FIRE TO REDUCE THE ATTITUDE 

* ANGLES TO ZERO. IF THE FLAG IS 1 THE CONTROL SYSTEM WILL REDUCE THE 

* NUTATION TO NUREF ONLY. 

IF ( (NU-MUREF ) .LE. 0.00001) THEN 
TX =0.0 
TZ = 0.00000 
ELSE 

IF ( FLAG. EQ.O) THEN 
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DNU = (NU-NUREF ) /2 . 0 + NUREF 
FLAG = 1 
END IF 

I F ( AY . GT . AON . AND . FLAG . EQ . 1 ) THEN 
FUME = TIME 
TZ = 2 . 0*F*R*FARC 
TX = 0.0 
ELSE 

TZ=0 . 0 
TX =0.0 
END IF 

IF ( MU. LE. DNU) THEM 

CANCEL = FTIME + ABS (3 . 14159/WN) 

FLAG = 2 
END IF 

* CALCULATION OF THE CORRECT THRUSTER COMBINATION REQUIRED TO REDUCE 

* THE ATTITUDE ANGLES TO A MINIMUM. 

I F (TIME . GE . CANCEL . AND . FLAG . EQ . 2 ) THEN 
DT = CANCEL + 2*FARC 
CONT = 4*F*R*FARC 

TCL = -CONT*STEP (CANCEL) + CONT*STEP (DT ) 

IF ( SIGMA. GT. 1.0) THEN 

BETA = ((SIGMA-1. 0)*3. 14159/SIGMA) -6.283 
ELSE 

BETA = (SIGMA-1. 0)*3. 14159/SIGMA 
END IF 

IF (3ETA. GT. 0.0. OR. BETA. LE. -1.5708) THEN 
DELTA = ABS (BETA) 

TZ = -TCL*COS (DELTA) 

TX = -TCL*S IN (DELTA) 

ELSEIF (BETA. GT. -1.5708. OR. BETA. LE. -3. 14159) THEN 
DELTA = 3.14159- ABS (BETA) 

TZ = TCL*COS (DELTA) 

TX = -TCL*S IN (DELTA) 

ELSEIF (BETA. GT. -3. 14159. OR. BETA. LE. -4. 7123) THEN 
DELTA = ABS (BETA) - 3.14159 
TZ = TCL*COS (DELTA) 

TX = tcl*sin(delta) 

ELSEIF (BETA.GT. -4. 7123. OR. BETA. LE. -6.283) THEN 
DELTA = 6.283 - ABS (BETA) 

TZ = -TCL*COS (DELTA) 

TX = TCL*S IN (DELTA) 

END IF 
END IF 
15 END IF 

CONTRL FINTIM=20 . 0 , DELT=0.01 

SAVE .01, NU, PHI, PSI, WX, WZ , WY, SIGMA, LAMDA , WN , AY, AON , . . . 

H, HS, HT, PHID, PSID, THD , TH, TX, TY, TZ, WT 
GRAPH (G1,DE=TEK618) TIME(UN=SEC) PHI (UN=RAD ,LO=-3 . 0E-03 , SC=1 . 5E-03) . . 

PSI (UN=RAD ,LO=-3 . 0E-03,SC=1.5E-03) 

GRAPH (G3,DE=TEK618) TIME NU(UN=RADS) 

GRAPH (G4,DE=TEK618) TIME TX TZ (L0=-32E-03 , SC= . 0165 ) 

LABEL (G1,DE=TEK618) ROLL AND YAW ANGLES; WY = 20RPM 

LABEL (G2,DE=TEK618) YAW ANGLE; SIGMA=0.3635 & 1.479; WY=-6.283 

LABEL (G3,G9,DE=TEK613) NUTATION ANGLE; WY = 20 RPM 

LABEL (G4,DE=TEK618) FIRING DIAGRAM WY = 20 RPM 

END 

STOP 
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